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Abstract 

The symmetry of exchange interaction of charge carriers in semiconductor 
nanostructures (quantum wells and quantum dots) is analysed. It is shown 
that the exchange Hamiltonian of two particles belonging to the same energy 
band can be universally expressed via pseudospin operators of the particles. 
The relative strength of the anisotropic exchange interaction is shown to be 
independent of the binding energy and the isotropic exchange constant. 
The reduced symmetry of semiconductor nanostructures suggests that the exchange in- 
teraction of charge carriers in such structures is not necessarily described by the isotropic 
(Heisenberg) spin Hamiltonian. In particular, the exchange interaction of electrons and 
holes in quantum wells and quantum dots is known to be extremely anisotropic, giving 
rise to a fine structure of nanostructure excitons 0. It has been shown p[ that the ex- 
change interaction of conduction-band electrons is also anisotropic if the structure lacks 
inversion symmetry. The main term of the anisotropic exchange Hamiltonian in this case 
has the Dzyaloshinskii-Moriya form [Q. The electron-electron anisotropic exchange was 
subsequently widely discussed in relation to the quantum computing problem It has 

been recently detected experimentally via its contribution to the spin relaxation of donor- 
bound electrons in GaAs [0, where it has been shown to put the upper boundary for the 



electron spin lifetime at donor concentrations around 10^^ cm~^ . However, it remained so 
far unclear whether or not the anisotropic spin Hamiltonian suggested in Ref. is universal 
for all types of charge carriers, e.g. for two-dimensional holes. The issue of the dependence 
of the anisotropic exchange constant on the parameters of the localizing potential is also 
very sensitive, especially for the discussion on feasibility of quantum computation with solid- 
state spin systems The constant was so far calculated using the Heitler-London method 
PJ^, which is known to give incorrect asymptotic expression for the isotropic exchange in- 
tegral HJ^. Gor'kov and Krotkov ||lO|, using the median-plane method [§, have recently 



obtained a correct asymptotic formula for the anisotropic exchange constant in a specific 
case of hydrogen-like centers in zinc-blende semiconductors, different from that calculated 
earlier by the Heitler-London method 0. However their approach is not always applicable 
to coupled quantum dots, where the distance between quantum dots can be comparable to 
the quantum-dot size. 

The collection of unsolved problems and blanc spaces in the existing knowledge on the 
anisotropic exchange in semiconductor structures, given above, demonstrates the evident 
demand for a consistent theoretical analysis of the issue, based on a general approach. In 
this paper, we consider exchange interaction of two identical charge carriers localized in any 
symmetric double-well potential in a two-dimensional semiconductor structure. Using the 
pseudospin formalism allows to obtain a universal spin Hamiltonian describing this class of 
systems. 

Let us consider the exchange interaction of two identical charge carriers (electrons or 
holes), localized in two centrosymmetric potential hollows (further referred to as quantum 
dots, QDs) in a quasi-two-dimensional semiconductor structure (quantum well, QW). The 



QDs may be, for example, self-organized QDs otherwise, they can be induced by electro 



static potential of nanometer-sized gates or impurity centers |jT^ . The distance between 
centers of the QDs will further be denoted as R12. In quasi-two-dimensional structures, the 
4-fold degeneration of the valence band, typical of cubic semiconductors, is lifted. The states 
at extremum points of two-dimensional subbands in absence of magnetic fields retain only 



the Kramers two-fold degeneration. Their wave functions can be written as \l/(r)Mi,(r), where 
is an envelope function, and Ui,{r) is a Bloch amplitude, u = ±1/2. The Bloch amplitudes 



Uu{r) are transformed into each other by the operator of time reversal [|I4 



M_i/2(r) = -i(3-yM+i/2(r) (1) 

This property allows to associate the Kramers index u with an eigenvalue of a projection 
of a pseudospin operator j (j = 1/2) on some (generally, fictitious) axis. The choice of basis 
functions for j is not unambiguous. It is limited only by the condition given by Eq.(l). In 
particular, for heavy holes with the projection of the angular momentum on the structure axis 
Z, equal to = ±3/2, it is convenient to choose the functions as |jl|: |j, ±1/2) = | J^, —3/2) 
and |j, — 1/2) = |J2,±3/2). This choice allows to avoid phase multipliers which would 
otherwise appear at wave functions in the pseudospin representation. For conduction-band 
electrons, the pseudospin coincides with the electron spin s. Linear transformations of 
pseudospin wave functions determined in the basis ^^-1/2} are equivalent to rotations 



of usual spinor functions [14 



'"+1/2 ~ exp(27/2) M+1/2 exp(ia/2) cos(/9/2) ± exp(— ia/2) sin(/3/2) 

■"-1/2 = exp(— 'j7/2) — ^+1/2 exp(ia/2) sin(/3/2) ± ■u_i/2 exp(— 'ia/2) cos(/?/2) (2) 

where a, (3, and 7 are analogs of Euler angles. Following the analogy, one can introduce 
the total pseudospin I = ji ± j2- Indeed, the Gilbert space of two-pseudospin wave functions 
Aij.uUn{ri)ui,{r2) breaks into two subspaces invariant with respect to the simultaneous trans- 
formation of both pseudospins along Eq.(2) with the same a, f3, and 7. The basis functions 
of these subspaces, ^0 = (M+i/2(ri)'U_i/2(r2) - 'U-i/2(ri)M+i/2(r2))/V2 and ^im, equal to 
(M+i/2(ri)M_i/2(r2) ±M+i/2(ri)M_i/2(r2))/V2 (M = 0) or M±i/2(ri)M±i/2(r2) (M = ±1), are, 
obviously, eigenfunctions of the operators P and J^. 

The general form of the one-particle Hamiltonian of the two dimensional charge carrier 
in the pseudospin representation is 

^i = |^P±V^(r)±h(k)-j (3) 



where the "spin-orbit field" h(k) is a vector in the pseudospin space |T^. h(k) is an 
odd function of the components of the wave vector. It is not equal to zero if the structure 
lacks inversion symmetry (which is very typical for nanostructures). This is the case when 
either the crystal unit cell lacks inversion symmetry (bulk inversion asymmetry, BIA [p^), 
or the QW is asymmetric (structure inversion asymmetry, SI A [|I7|). The components of 
h(k) may be, or may not be, associated with certain Cartesian axes in the real space. In 



the two-dimensional case, h(k) is dominated by linear in k terms |1T7|JT8| , [19 



hrj = Arj^k^ (4) 

where the matrix A is defined by the structure symmetry. 

The problem we are going to solve is finding the fine structure of the ground state of the 
two-particle Hamiltonian: 

^ = ^0 + Hso (5) 



where 



Ho = T^kj + —kl + Vi(ri) + V2{r2) + Uu{\ri - r^l) (6) 
2m 2m 



^so = h(ki)-ji+h(k2)-j2 (7) 

and f/i2(|i"i — 1*2 1) is the operator of the Coulomb interaction between the two particles. 
Before tackling the effects of spin-orbit interaction in the form of Eq.(7) on the exchange 
interaction, we should reconsider the ground-state structure of the Hamiltonian Hq. It is 
indeed well-known for electrons whose one-particle wave functions are \l/(r)(^^, where the 
spinor is not a function of coordinates. To the contrary, the Bloch amplitude u^, does 
depend on coordinates, and, moreover, it may contain spinors with both fi = +1/2 and 
/i = —1/2. The exciton (an electron-hole pair) in a QD is a good example demonstrating that 
the exchange interaction of charge carriers may have a very different symmetry as compared 
to that of free electrons. The QD exciton fine structure UTI] consists of two doublets, being 



thus quite different from the fine structure of a pair of vacuum electrons, i.e. the well-known 
singlet-triplet structure associated with the Heisenberg exchange. 

In order to analyze the fine structure of Hq for two particles belonging to the same 
subband, we first note that their behavior should be identical to that of bare electrons in all 
aspects but the Coulomb interaction. Indeed, although Bloch amplitudes are functions of 
coordinates within the unit cell, the one-particle operators of the kinetic energy and of the 
potential energy in Eq.(6) in the effective-mass approximation act upon envelope function, 
not Bloch amplitudes. Therefore, with respect to these one-particle operators, the Bloch 
amplitudes are just equivalent to spinors. To the contrary, calculating the Coulomb energy 
assumes taking integrals over the unit sell also. It is due this fact that the symmetry of 
Bloch amplitudes of holes and electrons has an impact on the fine structure in the exciton 

EJH- 

The fermionic wave functions of the two charge carriers can be written in the following 
form, similar to that of bare electrons: 

^o(ri, r2) = [$o(ri, rs) + $o(r2, ri)] 
^iM(ri, ra) = [$i(ri, rs) - $i(r2, ri)] (8) 

where $o(ri,r2) and $i(ri,r2) are two-particle envelope functions defined so that each 
particle is most likely to be found near its "home" center, while $o(r2,ri) and <l'i(r2,ri) 
correspond to interchanged particle positions. To determine the structure of respective 
energy levels, we should recall a property of the Bloch amplitudes u^, which results from 
their symmetry with respect to time reversal, and is an equivalent formulation of the Kramers 
theorem. As it follows from the Kramers theorem, the states symmetric with respect to time 
reversal remain degenerated unless magnetic field is applied. Mathematically, this means 
that matrix elements of any function of coordinates (not containing derivatives or spin 
operators) between u^i/2{r) and u^i/2{r) are zero, while diagonal matrix elements are equal 
to each other: 
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J n+i/2(r)n*_i/2(r)/(r)ciV = J u^„2{v)u\y^{v)f{v)dh = 

I «+i/2(r)<i/2(r)/(r)t/V = J u^„^{v)uU/2{^)f{v)dh (9) 

where the integrals are taken over the unit celL 
Using Eq.(9), one can easily find that 

(^o(ri, rs)! t/i2(|ri - ral) |^iM(ri, ra)) = 
(*o(ri,r2)|f/i2(|ri-r2|)|v&o(ri,r2)) ^ (10) 
(^iAf(ri,r2)| f/i2(|ri - r2|) |^iM'(ri5 ^2)) = const ■ 5mm' 

Thus, the Coulomb interaction retains the singlet-triplet structure of the ground state 
of two identical charge carriers. Exactly like in the case of two bare electrons, two-particle 
states with the same total pseudospin / are degenerated. Consequently, the Hamiltonian of 
the exchange interaction in terms of pseudospin operators takes the Heisenberg form: 

^5 = -2A(ji-j2 + l) (11) 

where A is a constant to be determined for each specific case. 

Now we can consider the effect of the spin-orbit terms given by the Eq.(7) on the exchange 
interaction. In the following, we will choose the axis X along the straight line connecting 
the localization centers (QDs). To handle the spin-orbit terms, we make use of a unitary 
transformation proposed by Levitov and Rashba who used it to eliminate spin-orbit 
terms in the one- dimensional case. The matrix T defined as 



2m 

>, — : 

h 



2 ^ AaxihaXl + 32aX2) 



(12) 



T = exp 

transforms the Hamiltonian Eq.(5) into the form: 

THT-^ = H' = H'q + H'so (13) 

where 
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and 

= E Jla + (15) 

a 

where j'l = TjiT~\ j'^ = T^T''^- 

The Hamiltonian H' does not contain spin-orbit terms and therefore results in the ex- 
change interaction in the form of Eq.(ll): 

H's = -2A(j;-j'2 + 1) (16) 

Due to the axial symmetry of the system, the matrix elements of kiy and k2y, calculated 
on the ground-state eigenfunctions of Hq (they can be obtained from Eq.(8) by the transfor- 
mation with the matrix T, which does not affect their dependence on y), are exactly equal 
to zero. The same is true for all the odd powers of kiy and k2y. Therefore, H'so does not 
contribute into the exchange interaction. 

Finally, to obtain the exchange Hamiltonian in the non-transformed basis, one should 
substitute the expressions for j'^ and j2 into Eq.(16). Since the transformation T is a rotation 



through the angle JY, ^vx^vxX around the vector Aj,x in the pseudospin space, Hs in 



the non-transformed basis is not unambiguously defined: it depends on the coordinates xi 
and X2 at which we take the spin operators ji and j2. A natural choice is to define them at 
the centers of corresponding QDs; for instance, this definition allows to write the Zeeman 
interaction in the usual form, Hz = fJ'BQa/sjaB/s, where B is the magnetic field, /x^ is the 
Bohr magneton, and gaf3 is a symmetric tensor g-factor p2| whose principal directions do 
not depend on the envelope wave function of the localized particle. This way, we come to 
the expression for Hs obtained in Ref. 0: 

^5 = -2A(l+jrj2COS7 + (d-ji)(d-j2)(l-cos7) + d-[ji X j2]sin7) (17) 



where 7 = ^ JJ2 ^uxA^xRu, and d is a unit vector in the pseudospin space, defined 



u 



so that = Ayxl . Yj ^vx-^vx- The first anisotropic term has the form of pseudodipole 



interaction |2^, and the second one, of the Dzyaloshinskii-Moriya interaction ||^. At small 
7, the Dzyaloshinskii-Moriya interaction dominates. 

The Eq.(17) demonstrates a remarkable universality of the exchange interaction in two- 
dimensional semiconductor nanostructures: this form of the Hamiltonian holds for both 
electrons and holes, for any type of centrosymmetric localizing potentials. The Eq.(17) is 
valid for identical as well as for different QDs. Moreover, the angle 7 characterizing the 
relative strength of the anisotropic exchange depends only on the distance between the QDs 
and the orientation of the pair of QDs with respect to the crystal axes. It is not sensitive to 
binding energies of the charge carriers in the QDs and to the value of the isotropic exchange 
constant A. 

The value of 7 can be now easily calculated for those structures where the components 
of the matrix A are known. 



In [100] oriented GaAs quantum wells the dominating BIA terms are Ayy = —A^ 



(kl), Axy = Ayx = where a ~ 0.65 [E^ (here coordinates x and y are taken along 



m 

the cubic crystal axes). This gives 7 = (A;^) i?i2. For example, in a 5nm-wide QW 

considered in Ref. 7 = (3 ■ 10^cm~^)i?i2, which gives for R = 24nm 7 = 0.72, i.e. more 
than an order of magnitude greater than the value of 0.02 obtained in Ref. at the same 
values of parameters by the Heitler-London method. 

For Rashba terms, Ayy = A^x = 0, A^y = —Ay^ = a. In a single-side modulation-doped 
n-type Si/SiGe quantum well, the constant a of 1.1-10"^^ eV-cm was measured This 
gives 7 = ^Ri2 ^ (6.7 ■ W^cm-^)Ri2. 

Bulk inversion asymmetry terms for holes in zinc-blende semiconductors include both 



cubic and linear in k terms |^5[. The cubic term H^y has the same symmetry as the 
Dresselhaus term for electrons, with the constant — ^^y= — ; where m is the conduction-band 

mw 2mEg 

electron mass, and ay ~ 0.1 for GaAs. The linear term is given by the expression: 

//ly = -lx(k.f2) (18) 

where = (^J^ — Jy^ + {Jl — Jy^ Jz (other components of f2 are obtained by cyclic 



interchange of indices), and k ^ 10~^°eV^ ■ cm. Taking matrix elements of }i\y and H-iy 
within pairs of the states with = ±1/2 (hght holes) and = ±3/2 (heavy holes), and 
going to the pseudospin notation, we obtain for a [100] QW: 

h,(k) = - [2^^= (kl) - 2V^X I {k^j, - kyjy) 

y ■mJ2mEg J 
h/j(k) = -iVSx^k^j^ - kyjy) (19) 

where h/(k) and h/i(k) are spin-orbit fields (see Eq.(3)) for light and heavy holes re- 
spectively. Consequently, 7; = ( (k^) - ^^x) Ru for light holes, and = 
sVsrrih ^ji^^ fQj. heavy holes, where mi and rrih are effective masses of the light and heavy 
hole respectively, corresponding to their motion along the QW plane. For example, in a 
lOnm-wide GaAs QW, 7^ ^ 7/, ^ (2 ■ l0^cm-^)Ri2. 

The symmetry of exchange interactions has been discussed in relation to feasibility of 
quantum computation with spins of localized electrons in semiconductor nanostructures 
0,^]^. A necessary (but not sufficient; see Ref. [^) condition for practical quantum com- 
puting to become possible is that the error probability per quantum gate be less than a 
certain value (of the order of 10"^) |^^. As shown in Ref. 0, there exists a way of perform- 



ing exchange-mediated quantum gates that allows to avoid errors caused by the anisotropy, 
provided 7 remains constant when the isotropic exchange constant is changed. The above 
consideration shows that this is indeed the case as long as spin-orbit constants do not 
change. They may change, however, because application of electric fields to the structure 
can either alter {kl) or bring about SIA (Rashba) terms. Since typical values of 7 are in 
reality much greater than estimated in Ref. 0, the uncontrollable effect of anisotropy on 
exchange-mediated quantum gates can not be so easily discarded. 

In conclusion, the exchange interaction of charge carriers (electrons or holes) localized in 
two-dimensional semiconductor structures is shown to be described by an universal Hamil- 
tonian in terms of carriers' pseudospins. It has the Heisenberg form unless spin-orbit terms, 
linear in the carrier wave vector, are present in the total Hamiltonian of the system. In this 
latter case, anisotropic contributions having both Dzyaloshinskii-Moriya and pseudodipole 



form arise. The "rotation angle" 7, characterizing the relative strength of the anisotropic 
exchange, linearly depends on the distance between the localization centers and does not 
depend on binding energies of the carriers. 
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